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ABSTRACT 


iCrl-d 


The  flow  field  in  a  two-dimensional  mixing  layer,  highly  disturbed  by 
a  sinusoidally  oscillating  flap,  was  mapped.  Two  velocity  components  were 

measured  throughout  the  flow  field  using  a  rake  of  X-wire  probes. 

/ 

Streaklines  were  calculated  from  the  phase^-locked  measured  data,  and  were 

compared  to  smoke  pictures,  creating  a  link  between  flow  visualization  and 

/ 

experimental  results.  The  phase^locked  vorticity  field  and  the  Reynolds 
stresses  were  calculated,  from  the  experimental  data. 
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It  was  found  that  fluctuations,  locked  in  phase  with  the  disturbance 
frequency,  were  not  only  responsible  for  the  fast  initial  growth  of  the 
mixing  layer  but  also  for  its  shrinkage  further  downstream  ^^the  occurance 
of  regions  I  and  II  in  the  parlance  of  Os ter  and  Wygnanski).  ''The  resump¬ 
tion  of  the  growth  of  the  mixing  layer  further  downstream  (region  III)  Is 
no  longer  controlled  by  phase-blocked  oscillations.  Although  these  oscilla¬ 
tions  still  exist  in  region  III  they  lose  energy  with  increasing  distance 
from  the  splitter  plate.  The  first  subharmonic  of  the  Imposed  frequency 
was  not  detected  and  the  pairing  of  vortices  in  regions  II  or  111"*  was  not 
observed. 

A  comparison  was  made  between  the  experimental  results  and  calcula¬ 
tions  based  on  linear,  inviscid  stability  theory  for  a  slowly  diverging 
plane  mixing  layer,  Good  agreement  was  observed  between  the  cross  flow  dis¬ 
tributions  of  any  predicted  quantity  and  a  measured  one,  however,  the  theo¬ 
retically  predicted  rates  of  amplification  with  distance  from  the  splitter 
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plate  were  higher  than  the  observed  ones.  The  linear  theory  Is,  therefore, 
Halted  In  Its  ability  to  model  the  entire  behavior  of  the  disturbed  mixing 
layer.  -  c  '  , 
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CHAPTER  1 


INTRODUCTION 


The  turbulent  mixing  layer  between  two  uniform  coflowing  streams  hav¬ 
ing  a  different  velocity  (Fig.  1.1)  has  been  the  subject  of  many  investi¬ 
gations  over  the  past  few  decades.  The  vast  interest  in  this  simple  confi¬ 
guration  is  undoubtably  trigerred  by  the  numerous  natural  phenomena  with 
which  it  can  be  associated  and  its  wide  range  of  applications. 

Two  approaches  are  currently  used  to  study  turbulent  shear  flows  of 
this  kind.  The  classical  approach  views  turbulence  as  a  random  process 
which  can  only  be  adequately  described  statistically.  The  other  approach 
looks  for  quasi-deterministic  vortex  structures  which  are  considered  to  be 
the  main  building  blocks  of  such  flows.  The  adoption  of  the  latter  appro¬ 
ach  in  the  investigation  of  the  two  dimensional  mixing  layer  was  stimulated 
by  flow  visualization  experiments  of  Brown  &  Roshko  (1971,1974)  who  disco¬ 
vered  the  presence  of  large  coherent  structures  in  the  turbulent  mixing 
layer  between  two  inert  gases.  Winant  &  Browand  (1974)  observed  a  merger 
between  pairs  of  adjacent  eddy  structures  to  which  they  attributed  the  la¬ 
teral  growth  of  the  turbulent  mixing  layer.  This  "pairing  process"  occurs 
randomly  in  space  and  time  (Roshko  1976)  generating,  on  the  average,  a  smo¬ 
oth,  linearlly  growing  mixing  layer.  It  is  now  generally  accepted  that  the 
large  eddies  play  a  major  role  in  the  growth  of  the  mixing  layer  by  influ¬ 
encing  the  entrainment  process;  by  amalgamating;  and  by  a  mutual  destruc- 


tion  referred  to  as  "tearing"  (Dimotakis  &  Brown  1976). 
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Reservations  regarding  the  two-dimensionality  of  the  large  eddies  and 
their  indefinite  persistence  may,  however,  still  exist.  Chandrsuda  et  al. 
(1978)  suggest,  that  the  two-dimensional  eddies  represent  a  relic  of  the 
transition  process,  while  the  fully  developed  turbulent  flow  is  highly 
three-dimensional.  They  argued  that  natural  transition  does  not  occur  whe¬ 
never  the  level  of  the  free  stream  turbulence  is  high  and  therefore,  large 
coherent  structures  can  only  appear  under  ideal  laboratory  conditions. th 
Pui  &  Garthshore  (1979)  associated  the  periodic  behaviour  of  the  large  ed¬ 
dies  with  vibrations  of  the  apparatus.  The  structures  observed  by  them 
were  not  even  moderately  two-dimensional.  Wygnanski  et  al.  (1979)  intro¬ 
duced  strong  external  disturbances  into  each  of  the  streams  (u'/U  ■  3*)  and 
found  quasi  two-dimensional  large  eddy  structures  forming  irrespectively  of 
the  high  turbulence  level  existing  in  the  free  streams.  Breidenthal  (1980) 
disturbed  the  flow  in  the  spanwise  direction  by  using  a  specially  designed 
splitter  plate  and  found  that  the  shear  layer  quickly  forgets  the  those 
perturbations  and  forms  the  characteristic  two-dimensional  vortex  struc¬ 
tures.  Nevertheless,  streamwise  vortices  which  appear  to  scale  with  the 
spanwise  vortices  persist  under  some  conditions.  (Roshko,  private  communi¬ 
cation) 

Observations  showing  the  influence  of  periodic,  two-dimensional  exci¬ 
tation  on  the  mixing  layer  were  reported  by  Oster  et  al.  (1978)  and  later 
by  Oster  6  Wygnanski  (1982)  .  The  rate  of  growth  of  the  shear  layer  was 
significantly  altered  by  the  excitation  even  when  the  amplitude  of  the 
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latter  was  small.  Ho  &  Huang  (1980,  1982)  excited  the  laminar  mixing  layer 
at  various  frequencies  and  showed  that  adjacent  vortices  can  either  pair 


with  one  another  or  "collective  interact”.  Vortex  amalgamations  were  ar¬ 
tificially  stimulated  by  these  investigators  who  prompted  the  coalescence 
of  any  small  number,  N,  of  vortices  by  disturbing  the  flow  with  the  Nik* 
subharmonic  of  the  initially  most  amplified  frequency. 


Statistical  theories  in  which  the  turbulence  was  considered  to  be  a 
random  process  were  always  impeded  by  the  fact  that  the  number  of  unknowns 
exceeded  the  number  of  equations.  One  had  to  resort  to  models  (e.g.  eddy 
viscosity,  mixing  length,  turbulent  rate  equation,  a  relation  between  Rey¬ 
nolds  stress  and  turbulent  Intensity  etc.)  to  provide  the  additional  equa¬ 
tion  required  for  the  solution.  The  discovery  of  large  coherent  structures 
suggest  new  ways  of  averaging  and  entirely  different  approaches  to  model¬ 
ling.  The  mixing  layer  is  sometimes  replaced  by  an  array  of  vortices 
evolving  in  time  (Corcos  &  Sherman  (1976),  Acton  (1976),  Ashurst  (1976), 
Patnaik,  Sherman  &  Corcos  (1976),  Corcos  &  Sherman  (1984),  Corcos  &  Lin 
(1984)  and  Lin  6  Corcos  (1984),  to  name  a  few). 


With  the  Increased  availability  of  large  computers  numerical  solutions 
of  the  equations  of  motion  became  feasable  using  an  ever  decreasing  number 
of  constrains.  Large  eddy  simulations  were  applied,  among  others,  by  Riley 
&  Metcalfe  (1980).  They  found  that  the  presence  of  the  subharmonic  fre¬ 
quency  which  is  not  in  phase  with  the  most  unstable  linear  eigenmode, 
causes  the  vortices  to  pair.  Riley  &  Metcalfe  suggested  a  plausible  expla¬ 
nation  for  some  of  the  observations  of  Oster  et  al.  (1978)  in  particular 
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the  discontinuity  in  the  rate  of  growth  of  the  mixing  layer  which  was  ac¬ 
companied  by  the  appearance  of  negative  Reynolds  stresses.  They  explained 
the  effect  by  eliminating  the  subha nnonic  frequency  from  their  calculations 
which  resulted  in  the  reduction  of  the  turbulent  kinetic  energy  and  a  con¬ 
comitant  slow-down  in  the  lateral  rate  of  spread  of  the  flow.  It  was  sug¬ 
gested  that  the  elimination  of  the  subharmonic  inhibited  or  even  eliminated 
the  pairing  process. 

Inoue  (1983)  simulated  vortex  pairing  in  turbulent  mixing  layers  by  a 
discrete  vortex  method.  He  found  two  stages  in  the  developement  of  the 
mixing  layer:  a  stage  of  clustering  where  numerous  vortices  interact  with 
one  another  and  form  discrete  lumps  followed  by  a  pairing  stage  in  which 
the  lumps  amalgamated.  Inoue  compared  his  calculations  with  experiment  and 
found  reasonable  qualitative  agreement  with  available  data. 

The  evolution  of  the  large  eddies  can  also  be  predicted  from  hydrody¬ 
namic  stability  theory.  In  laminar  flows,  where  the  mean  velocity  profile 
represents  the  solution  of  the  equations  of  motion,  one  may  determine  the 
character  of  the  disturbances  resulting  from  the  instability  of  the  flow 
field.  Such  an  analysis  was  done  by  Michalke  (1964,  1965)  who  predicted 
the  temporal  and  spatial  growth  of  perturbations  in  an  inviscid 
two- d omens ional  mixing  layer.  The  mean  velocity  profile  for  which  calcula¬ 
tions  were  performed  was  represented  by  a  hyperbolic-tangent  function. 
Monkevltz  and  Huerre  (1982)  compared  the  stability  of  the  "tanh”  profile 
with  the  profile  resulting  from  the  solution  of  the  boundary  layer  equa¬ 
tions  while  studying  the  influence  of  the  velocity  ratio  on  the  spatial 
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instability  of  mixing  layers.  They  found  that  the  maximum  amplification 
rate  of  the  perturbation  is  approximately  proportional  to  the  parameter 
^2~Ui)/(^2+^1  )  t*ie  shape  of  the  velocity  profile  was  less  important. 
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The  application  of  the  stability  theory  to  turbulent  flows  is  not 
straightforward  because  there  is  no  steady  velocity  field  in  existance  upon 
which  a  perturbation  may  be  superimposed.  One  tends  to  select  the  mean 
velocity  profile  as  being  representative  on  the  average  recognizing  the 
possibility  that  at  any  instant  in  time  the  actual  flow  field  may  be  quite 
different.  Nevertheless,  knowing  that  the  random  changes  in  the  mean  velo¬ 
city  occur  on  a  time  scale  which  is  short  in  comparison  with  the  period  as¬ 
sociated  with  the  large  coherent  structures  and  assuming  that  the  kinetic 
energy  exchange  between  largely  disparate  scales  is  negligible,  Gaster  Kit 
&  Wygnanski  (to  appear  in  J.F.M  )  applied  the  inviscid  linear  stability 
theory  with  a  large  measure  of  sucess.  They  modelled  the  development  of 
normal  modes  in  a  slowly  divergent  turbulent  mixing  layer  by  expanding  the 
inviscid  equations  in  multiple-scales  in  a  way  similar  to  the  one  used  by 
Bouthier  (1972)  and  by  Crighton  &  Gaster  (1976).  The  streamwise  velocity 
component  of  the  computed  eigenfunctions  agreed  well  with  experiment  but 
the  quantitative  comparison  of  the  growth  rates  was  poor. 


One  of  the  aims  of  the  present  work  was  to  extend  this  comparison  to 
the  lateral  velocity  component,  the  phase-locked  u'v'  product  and  the  span- 
wise  vorticity  distributions  and  to  determine  the  limitations  of  the  linear 


approach  by  applying  relatively  high-amplitude  disturbances  to  the  flow. 
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Furthermore,  most  of  the  observations  related  to  large  coherent  struc¬ 
tures  are  based  on  flow  visualization.  In  some  instances  the  method  of  vi¬ 
sualization  represents  local  conditions  (shadowgraph,  laser  induced  floure- 
sance,  chemical  reaction  etc.)  in  others  it  is  based  on  the  interpretation 
of  streaklines  (hydrogen  bubbles,  smoke  or  dye  filament  etc.).  The  often 
used  streakline  method  may  be  misleading  because  a  congregation  of  marked 
particles  in  one  location  in  the  flow  may  be  the  outcome  of  cumulative 
forces  acting  over  a  large  distance  and  may  not  nacessarily  represent  a 
local  event  of  dynamical  significance.  The  dynamical  significance  of  vor¬ 
tex  pairing  is,  for  the  time  being,  unknown.  We  therefore  tried  to  calcu¬ 
late  streaklines  from  phase-locked  velocity  measurements;  if  the  evaluated 
streaklines  agree  with  those  observed  by  marking  fluid  particles  then  the 
same  data  can  be  used  for  the  evaluation  of  phase-locked  vorticity,  u'v' 
product  and  any  other  quantity  of  siglnificance.  Therefore,  a  second  aim 
of  the  present  investigatin  was  to  develop  a  tool  which  does  not  rely  on 
qualitative  interpretation  of  marked  particles. 
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CHAPTER  2 


THE  EXPERIMENTAL  FACILITIES  AND  EQUIPMENT 


The  mixing  layer  facility  in  Tel-Aviv  University  consists  of  two  inde¬ 
pendent  cascade  blower  wind  tunnels  dischargigng  into  a  common  test  section 
(Fig.  2.1).  The  two  tunnels  which  are  a  mirror  image  of  one  another,  are 
separated  initially  by  a  splitter  plate  which  extends  upstream  through  the 
contraction  section  and  into  the  settling  chamber.  The  splitter  plate  ends 
20  cms.  downstream  of  the  contraction  allowing  the  two  streams  to  become 
parallel  before  the  initiation  of  the  mixing.  The  trailing  edge  of  the 
splitter  plate  was  milled  at  an  included  angle  of  3°. 

Each  tunnel  consists  of  an  air-supply  system,  a  diffuser,  a  settling 
chamber  and  a  nozzle.  The  air  is  supplied  to  each  tunnel  by  a  backward 
facing  step  blower  driven  by  a  separate  D.C.  motor.  One  of  the  motors  has 
a  power  rating  of  15  KW  and  is  manufactured  by  Elekt romotoren  Uerke  Karl- 
kaiser  (Typ  GV  132. 24S),  it  is  controlled  by  a  DFVLR  DMR  II  digital  con¬ 
troller.  The  other  motor  is  a  5  HP  "Kinematic"  D.C.  motor,  controlled  by 
a  General  Electric  "Thymotrol"  control  system,  equipped  with  a  tachometer 
generator.  The  air  speed  of  the  lower  stream  can  change  continously 
between  0-13.5  m/sec  and  in  the  upper  stream  between  0-9  m/sec.  The  two 
independent  blowers  enable  one  to  work  at  any  velocity  ratio  desired,  over 
a  wide  range  of  velocities.  The  blowers  were  vibration-isolated  from  the 
rest  of  the  structure  and  equipped  with  filters  at  their  inlet.  A  small 
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settling  chamber  and  a  gauze  are  situated  between  the  blower  and  the 
diffuser  in  order  to  provide  resistance  and  equalize  spatially  the  flow 
entering  the  diffuser. 

Each  diffuser  is  subdivided  into  four  smaller  channels  having  an  equi¬ 
valent  cone  angle,  that  does  not  exceed  5°.  A  deep  cell  honeycomb  and  3 
turbulence  damping  screens  are  installed  in  the  settling  chamber.  The  con¬ 
traction  ratio  of  each  nozzle  is  7.3:1  and  a  representative  turbulence 

i 

level  at  the  exits  was  measured  to  be  0.402  in  the  upper  tunnel  at  3  m/sec, 
and  0.272  in  the  lower  tunnel  at  5  m/sec. 

The  test  section  is  2000  mm  long,  500  mm  high  and  bQO  mm  wide.  The 
side  walls  which  are  made  of  transparent  plexiglass  sheets  are  removable 
for  photographic  purposes  (l.e.  avoiding  reflections).  The  top  and  bottom 

i 

J  walls  are  made  of  flexible  P.V.C.  sheets  mounted  on  screw  jacks,  enabling 

one  to  control  the  streamwise  pressure  gradient. 

U  The  traversing  mechanism  supporting  the  measuring  sensors  has  three 

degrees  of  freedom: 

a)  Translation  in  the  streamwise  direction  (X-axis),  performed  manually 
by  a  sliding  wheel-cart  with  a  resolution  of  1  mm. 

b)  Translation  in  the  vertical  direction  (Y-axis),  performed  by  a  lead 
screw  having  a  4  mm  pitch. 

c)  Rotation  around  the  Y-axis,  performed  by  a  shaft  passing  through  a 
hollow  column  of  a  symmetric  airfoil. 
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Both  the  lead  screw  and  the  shaft  are  activated  by  separate  Slo-Syn 
M063FDU6  stepping-motors,  having  an  angular  resolution  of  1.8°  per  step, 
equivalent  to  translational  resolution  of  0.02  mm.  The  stepping-motor  ro¬ 
tating  the  probes  about  the  Y  axis  is  connected  to  the  shaft  through  a  gear 
having  a  reduction  ratio  of  10:1,  therefore  increasing  the  angular  resolu¬ 
tion  up  to  0.18°. 

Both  stepping-motors  are  controlled  by  a  Rapac  Electronics  2-axes  PSM 
15 IB  controller,  equipped  with  manually  and  externally  triggered  modes. 
The  external  operation  mode  enables  us  to  computer-control  the  movements. 

A  thin  flap  pivoted  at  its  leading  edge  generated  the  required  pertur¬ 
bations.  The  flap  is  10  mm  and  0.5  mm  thick  and  it  spans  the  entire  test 
section  along  the  trailing  edge  of  the  splitter  plate.  The  gap  between  the 
splitter  plate  and  the  flap  was  sealed  by  an  adhesive  tape.  The  sinusoidal 
oscillations  were  provided  by  two  voice  coils  connected  in  parallel,  acti¬ 
vated  by  a  Krohn-Hyte  5100A  function  Generator,  insuring  full  control  over 
the  frequency  and  amplitude  of  the  fluctuations.  The  signal  was  amplified 
by  a  AO  watt  Dynaco  Stereo  80  power-amplif ier.  The  sine-wave  signal  was 
simutaneously  processed  by  a  home-made  frequency-dividing  circuit,  creating 
a  synchronized  square-wave  subharraonic  signal  used  as  a  t rigger- source  for 
the  data  acquisition  process.  A  rake  of  7  X-wires  operated  at  a  constant 
temperature  mode  was  used  throughout  the  investigation.  Each  wire  is  made 
of  Tungstan,  5  micron  in  diameter,  welded  to  a  pair  of  prongs  1  mm  apart, 
creating  an  angle  of  45°  with  the  axis  of  the  probe.  The  length  to  diame¬ 
ter  ratio  of  each  wire  is  1/d  *280  and  the  frequency  response, as  determined 


LrO  • 


•  V  V  ’/  V  /  /  /.  t  •  w  1 *  «  •  .  »  .r.  « 


by  a  square  wave  method,  was  approximately  10  KHz.  The  multi-channel  con¬ 
stant  temperature  anemometers  were  designed  and  constructed  in  the  elec¬ 
tronics  workshop  of  the  Faculty  of  Engineering  in  Tel  Aviv  University.  The 
over-heat  ratio  was  set  to  be  1.5. 
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CHAPTER  3 


THE  CALIBRATION  METHOD 


The  calibration  method  for  an  array  of  X-wlre  probes  described  below 
is  only  possible  with  the  aid  of  a  computer. 

It  is  assumed  that  the  output  voltage  of  each  wire  in  the  X-probe  is  a 
function  of  the  modulus  of  the  velocity  vector  Q  and  the  angle  this  vector 
makes  with  the  axis  of  the  probe  a  : 

(Q»o)  flnd  ^2  *■  E£  (Q,a). 

It  is  also  assumed  that  within  the  range  of  calibration,  the  opposite 
is  also  true:  for  each  pair  of  Q  and  a,  only  one  pair  of  voltages  is  pos¬ 
sible,  or  that  a  given  pair  of  voltages  and  E2>  determines  the  velocity 
uniquely. 

Q  -  Q(E1,E2)  and  *-  a(El,E2). 

The  functions  Q  and  a  are  assumed  to  be  forth-order,  two-dimensional 
polynomials  of  the  form: 

Q  -  AiE^A2E^A3eJe2+A4eJE^  .  .  .  +A15 

and  a  similar  expression  for  a«  In  order  to  determine  the  15  coefficients 
for  each  polynomial,  a  two-dimensional  Least  Square  Fit  (ISFP)  procedure 
was  employed.  Meaurements  of  the  mean  voltages  E^  and  E2  were  taken  at  7 
different  air  velocities  and  11  different  yaw  angles,  between  -27°<a<27°, 
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for  each  velocity.  Hence,  77  different  pairs  of  and  formed  the  basis 
for  the  LSFP  procedure. 

The  operation  of  the  calibration  procedure  was  almost  fully 
computer-controlled.  The  digital  output  channel  of  the  computer  was  used 
as  a  trigger-source  for  activating  the  stepping-motor  controller  rotating 
the  probe  holder  and  changing  its  yaw  angle.  Data  acquisition  started 
right  after  the  movement  was  completed.  The  whole  calibration  procedure, 
for  a  rake  of  7  X-wire  probes,  lasted  no  longer  than  5  minutes. 
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DATA  ACQUIS  ITION  AND  NUMERICAL  CALCULATIONS 


Data  acquisition  and  numerical  analysis  were  performed  on  a  PDP  11/60 
minicomputer  equipped  with  an  analog  to  digital  converter  which  is  a  part 
of  the  LPA11  system.  The  computer  has  a  total  memory  of  128K  words  and  is 
operated  by  a  RSX-11M  operating  system  (version  3.2).  The  computer  is  also 
equipped  with  two  Digital  RK07  disc  drivers  and  two  Kennedy,  Model  9300, 
magnetic  tape  drives.  The  tape  drives  were  used  to  record  in  real  time  the 
output  voltages  of  the  hot-wire  anemometers  and  amplifiers.  Since  no  digi¬ 
tal  processing  was  carried  out  during  data  acquisition,  the  experiment  was 
shortened  significantly. 


The  A/D  converter  resolution  is  2.44  millivolts  per  digit.  The  sam- 

|  pling  frequency  was  2.048  KHz  per  channel  and  since  all  16  A/D  channels 

1 

1 

'  were  occupied  (14  hot  wire  channels  +  disturbance  channel  +  subharmonic 

square-wave  channel),  the  sampling  frequency  of  the  computer  was  32.768 
KHz,  far  below  its  80  KHz  limitation. 


The  first  step  in  the  numerical  computations  was  to  translate  the  vol¬ 
tages,  recorded  on  the  magnetic  tape  during  data  acquisition,  to  the  two 
velocity  components,  using  the  calibration  coefficients.  This  procedure 
produced  the  instantenous  values  of  Q  and  a,  and  the  velocity  components  u 
and  v  were  calculated  from  the  following  relations: 


u  =  Qcos(a) 


and 


v  -  Q8in(a) 


4.1.  Phase-locked  averaged  velocity  components 

The  phase-locked  averaged  velocity  vectors  were  calculated  from  the 
sampled  velocity  vectors  in  order  to  analyze  the  periodic  flow  phenomena 
having  the  disturbance  frequency,  or  its  first  subharmonic. 


The  numerical  representation  of  the  phase-locked  averaged  signal  is: 

<u(  1)>“  Su(*n+(j-l)At+i*NPP);  Where  u  is  the  original 

i-i  u 

vector  containing  the  velocity  values  sampled; 


At  is  the  time  interval  between  data  points; 


NPP  is  the  number  of  data  points  per  period  (T)  of  the  frequency  analyzed 
(i.e.  NPP-  T/At); 


*q  is  anarbitrary  initial  phase; 

and  N  is  the  number  of  time-periods  of  the  analyzed  frequency  contained 


within  the  original  velocity  vector. 


The  determination  of  and  NPP  required  the  simultaneous  sampling  and 
recording  of  both  the  sine-wave  signal  activating  the  flap,  and  its  first 
subharmonic  signal,  produced  by  the  electronic  frequency-dividing  circuit, 
designed  to  produce  a  synchronized  square-wave  signal. 


The  sketch  bellow  illustrates  the  pahse-locked  averaging  method. 
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The  initial  phase  points  $q  corresponds  to  a  point  at  which  the  squa- 
rewave  crosses  zero  with  a  negative  slope;  NPP  was  calculated  as  the  mean 
number  of  data  points  between  two  ad.ja  ent  initial-phase  points  (trace 
[B]). 


The  result  of  an  averaging  process  of  a  signal  [C]  dominated  by  a  cer¬ 
tain  frequency  (a  fundamental  frequency)  which  is  phase-locked  with  its 
first  subharmonic,  is  a  signal  [ D ]  containing  two  consecutive  waves  of  the 
fundamental  superimposed  on  the  subha rmonic.  All  the  random  high  frequency 
influence  (turbulent  behaviour)  is  filterred  out.  By  averaging  with  res¬ 
pect  to  the  fundamental  only  one  obtains  trace  [E].  Trace  [E]  is  also  ob¬ 
tained  whenever  the  two  frequencies  (i.e.  the  fundamental  and  its  subhar- 
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monic)  are  not  locked  in  phase. 

The  length  of  each  sampled  velocity  vector  was  512  points  (equivalent 
to  a  period  of  0.25  sec),  and  100  such  vectors  were  acquired.  Thus,  the 
average  phase-locked  to  the  subharmonic  frequency  was  calculated  from  400 
events. 

Other  flow  parameters,  such  as  the  mean  velocity  profiles,  the  R.M.S. 
values,  the  Reynolds  stresses  and  the  spectra,  were  calculated  from  the  or¬ 
iginal  time  series. 

4.2.  Reynolds  stress 

Three  different  ways  of  averaging  the  -u7v7  product  were  employed: 

a)  The  Reynolds  stress  is  a  temporally  averaged  product  of  the  two 

time  dependent  velocity  vectors: 

t(X,Y)  -  -p(u'(X,Y,t)*v'(X,Y,t)) 

where: u'(X,Y,t)«u(X,Y, t)-U(X,Y) 

v'(X,Y,t)-v(X,Y,t)-V(X,Y) 

b)  The  phase-locked  ensemble- averaged  -<u,v,>  product,  producing  a 
value  containing  the  influence  of  all  frequencies  which  are  locked  in  the 
phase  to  the  excitation: 

<T>  -  -PZT^v^ 


c)  The  spectrum  of  the  Reynolds  stress.  In  order  to  obtain  both  the 
power  spectra  and  the  spectra  of  the  u'v'  product,  we  have  used  a  fast  Fu- 
rier  transform  operator  (FFT).  The  result  was  a  vector  containing  both  the 
real  and  Imaginary  coefficients  of  various  frequencies  for  each  velocity 
component. 

-£(u)  -  C^F)  where  j)-Ax(  jJ+iB^  j) 

£(v)  -  C2(F)  where  C2( j)-A2(  j)+iB2( J) 

where  A(j)  and  B(j)  are  the  sine  and  the  cosine  coefficients  of  the  j1*1 
frequency,  respectively.  The  sampling  frequency  (2.048  KHz),  and  the 

length  of  the  vector  containing  the  velocity  signal  (512  data  points),  de¬ 
fine  the  lowest  frequency  measured,  and  the  spectral  resolution  to  be  equal 
to  4  Hz,  and  the  highest  frequency  to  be  equal  to  1.024  KHz. 

The  power  spectra  were  calculated  by  multiplying  each  of  the  coeffi¬ 
cients  in  C^(f)  by  its  complex  conjugate  C1(F)  ; 

Su2(F)-C1(F)*cJ(F)  where  Su2( J)-Aj( J)2+BL( J)2 

The  Reynolds  stress  spectra  were  calculated  from  the  real  part  of  the 
cross-spectrum  of  the  two  velocity  signals,  i.e.  the  product  of  one  velo¬ 
city  component  with  the  complex  conjugate  of  the  other  velocity  component: 

Suv(F)  -  0.5  REAL[C1(F)*C*(F)] 
where  Suv(j)  -  -0.5[A1( j)*A2( j)+BL( j)*B2( j)) . 

4.3.  Energy  production 

The  overall  energy  production  term  defined  by  -u’v*^  can  ^e  divided 
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into  the  phase-locked  (or  coherent)  energy  production  which  is  responsible 
for  energy  transfer  from  the  mean  flow  to  the  phase-locked  fluctuations, 
and  the  random  (or  turbulent)  energy  production.  The  overall  energy  pro¬ 
duction  terms  were  computed  by  multiplying  the  mean  -u*v*  product  with  the 
local  mean  velocity  gradient,  and  the  phase-locked  energy  production  was 
calculated  applying  the  same  procedure  to  the  phase-locked  -<u’v’>  product. 
The  mean  velocity  profile  employed  was  a  modified  hyperbolic  tangent  pro¬ 
file,  fitted  to  the  experimental  data  (for  details  see  chapter  6). 

In  addition,  we  have  calculated  the  local  integrals  of  the  overall  and 
phase-locked  energy  production  over  the  cross-section  of  the  flow, 
i.e.  _£(-u~v'~)dY  and_/(<-u’v’>g^)dY,  respectively. 

4.4.  Vorticity  field 

The  phase-locked  averaged  spanwise  vorticity  field  was  computed  by 
numerical  differentiation  of  the  phase-locked  averaged  velocity  vectors. 

The  resulting  phase-locked  averaged  vorticity  signal  <A(X,Y,t)>  was 
used,  both  for  plotting  the  isodynes  (isovorticity  contours)  and  for  furth¬ 
er  mathematical  operations  (i.e.  Fourier  transformation  at  the  disturbance 
frequency  in  order  to  obtain  the  vorticity  eigenfunctions  and  compare  them 
with  the  results  of  the  stability  theory). 


•v.y  v.y.  -v 

f  V.  -■  /  - 


''.OV. 


4.5.  Streaklines 


A  streakline  connects,  by  definition,  all  fluid  particles  that  have 
passed  through  a  fixed  point  In  space  (a  source  point),  during  a  given  per¬ 
iod  of  time. 

In  order  to  establish  and  draw  a  streakline  from  the  velocity  vectors 
given  In  the  Eulerlan  representation,  one  must  transfer  it  to  the  Lagrangi- 
an  representation,  i.e.  calculate  pathlines  given  by: 

r  -  r(R, t) 

where  R  Is  the  position  vector  of  the  particle  at  t*0. 


<«! 


Providing  the  Jacobyan  of  the  matrix  of  the  given  system  of  equations 
has  a  non-zero  value,  the  following  relation  can  be  written: 

R-R(r,t) 

A  streakline  at  a  given  time  t*  can  be  directly  calculated  by  the  fol¬ 
lowing  expression: 

r  -  r[R(r0,t0),t*] 

where  Tq  is  the  position  vector  of  a  source  point; 

h 

t  is  a  parameter  defining  the  time  at  which  the  atreakline  is  to  be  plot¬ 
ted; 

,  * 
and  t0  is  the  time  variable  0  <  t^  <  t  . 

Another  way  of  representing  a  pathline,  in  a  two-dimensional  flow 


field,  is: 
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u-dX/dt  ;  v*dY/dt 

with  X(t*tg)»XQ  and  Y(t*tQ)“YQ.  Xq  and  Yq  are  the  space  coordinates  of  a 
source-point. 


We  restricted  our  discussion  to  a  row  of  source-points,  displaced 
vertically  at  the  first  measuring  X-location  (X“200  mm).  Thus  we  can  re¬ 
present  the  pathlines  as  : 

rl  =  rl<r0»t0»t) 

where  r^  is  the  position  vector  of  the  source-points; 
ty  is  the  lntial  time  parameter 
and  t  is  the  time  variable. 

Since  the  flow  field  considered  is  periodic,  with  a  period  T,  it  follows 
that: 

rl(r0»t0»t^  *  ri(  ro»to+nT,t)  n-i,2 . 

The  pathlines  were  calculated  by  numerical  integration  of  the 
phase-locked  averaged  velocity  vectors. 

Each  source  point  r^  produced  NPP  pathlines  (i.e.  the  number  of  data 
points  in  the  phase-locked  averaged  velocity  vector).  The  pathline  vectors 
were  identified  by  an  index  j  defined  as  j-l+ty/At  where  ty  is  the  time  the 
particle  passed  through  the  source  point  and  At  is  the  time  interval 
between  data  points. 

Since  At  is  determined  by  the  sampling  frequency  the  integration  led 
us  to  space  coordinates,  other  than  our  measuring  locations.  Interpolation 
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techniques  were  employed  in  order  to  establish  the  velocity  values  at  those 
points.  A  spline  interpolation  procedure  was  first  used  in  the  transverse 
direction  (Y-  axis),  and  a  linear  interpolation  method  was  then  applied  in 
the  streamwise  direction  (X-axis).  Since  the  period  T  corresponds  to  NPP 
data  points,  all  time  considerations  reduce  to  a  proper  indexing  of  the  po¬ 
ints  in  the  phase-locked  averaged  velocity  vectors. 


The  function  r^  describes  the  position  of  all  particles  that  had 
passed  through  a  source  point.  Writing  : 

*\ 

ri“ri<ro*to,t  ' 

where  r^  and  t  are  now  the  parameters  of  the  source  position  and  the  plot- 

* 

ting  time,  respectively,  and  tQ  is  the  variable,  (  0<  ty<t  ),  we  thus  re¬ 
present  the  desired  streaklines.  The  periodic  flow  field  demands  that: 

^(ro.tQ.t  )  ■  +nT)  n*l,2 . 

Thus,  the  computation  of  streaklines  is  reduced  to  a  selection  procedure  of 
data  points  from  the  pathlines. 

If,  for  example,  the  streakline  at  the  time  corresponding  to  the  i**1 
point  of  the  phase-locked  averaged  vector  is  to  be  plotted  (i.e.  t*»i*At), 

the  following  points  are  to  be  picked  up  from  the  pathline  vector  indexed 

J-l: 

ith,  (NPP+i)th,  (2*NPP+  i)th,.  .  .  ,(n*NPP+i)th, 
from  the  pathline  indexed  j*2  the  following  points  are  to  be  picked  up: 

(i-l)th,(NPP+i-l)th,(2*NPP+i-l)th,  .  .  .  ,(n*NPPfi-l)th  , 

from  the  pathline  Indexed  j*j  the  following  points  are  to  be  picked  up: 
(i-JH)th,(NPm-j+l)th,(2*NPP+i-j+l)th,  .  .  .  ,(n*NPP+i-j+l)th  , 
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CHAPTER  5 


THE  DISTURBED  MIXING  LAYER 


5.1.  Mean  flow  parameters 

"By  introducing  a  periodic  disturbance  to  the  mixing  layer,  one  can 
change  dramatically  the  mean  velocity,  as  well  as  the  turbulent  properties 
of  the  flow"  (Oster  and  Wygnanski  1982). 

The  velocity  ratio  between  the  two  streams  in  the  present  experiment 
was  maintained  at  0.6,  the  velocity  of  the  slower  stream  was  U^-6  m/sec  and 
the  velocity  of  the  faster  stream  was  U2“10  «/sec.  The  flow  was  disturbed 
by  a  flap  oscillating  at  the  frequency  of  f-44.5  Hz,  and  an  amplitude  of 
1.5  mm.  At  each  cross  section  the  mean  velocity  profile  was  measured,  and 
from  it  the  mixing  layer  width,  b,  was  estimated  as  follows: 

b“(Y0.95”Y0.10) 

where  Yq  and  Yq  are  the  lateral  coordinates  at  which 

(U-U^ )/(U2*U^ )"0.95  ans  0.10  respectively.  The  momentum  thickness  9: 

m 

^-(U^)'2  J(U2-U)(U-U1)dY 

-it 

was  also  computed.  Both  quantities  are  plotted  against  the  streamwise 
coordinate  X  in  Fig.  5.1.1.  The  dimensional  streamwise  coordinate  X  is 
presented  to  avoid  ambiguity  resulting  from  the  dependence  of  6  on  X. 

YG.95’  Y0  50  dn<*  Yu  10  are  plotted  against  X  in  Fig.  5.1.2.  A  dense  grid 
of  data  points  (80  in  the  streamwise  direction)  was  necessary  to  achieve 
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sufficient  resolution  for  the  calculation  of  streaklines.  Unlike  the  un¬ 
disturbed  mixing  layer,  the  excited  mixing  layer  does  not  spread  linearly 
with  X  .  Three  distinct  regions  can  be  recognized  in  this  case  : 

Region  I,  (X<5b0  mm),  in  which  the  momentum  thickness  grows  with  incre¬ 
asing  distance  from  the  splitter  plate.  Part  of  region  1  (X<500  mm)  is 
characterized  by  a  linear  growth  of  the  momentum  thickness  with  a  slope  of 
dO/dX-0.0265  . 

Region  II , (560<X<840  mm),  in  which  the  width  of  the  mixing  layer  actual¬ 
ly  decreases. 

Region  III,(X>840  mm),  in  which  the  mixing  layer  resumes  its  linear 
growth,  but  at  a  slower  rate  (d©/dX*0.0071). 

Two  transition  subregions,  discussed  later,  separate  these  clearly  de¬ 
fined  regions.  This  behaviour  is  characteristic  of  a  strongly  disturbed 
mixing  layer  (Os ter  and  Wygnanski  1982). 


The  mean  velocity  profiles  scaled  by  the  momentum  thickness  are  simi¬ 
lar  in  both  regions  I  and  III.  The  actual  shape  of  the  profiles  changes 
from  region  to  region,  particularly  on  the  low-speed  side  of  the  flow  (Fig. 

5.1.3)  where  (i  ")  T>(^  ,  ttt'  The  difference  in  the  shape 

9  dY' region  I  0  dY' region  III  r 

of  the  mean  velocity  profiles  can  not  be  attributed  to  an  error  and  must 


stem  from  the  dynamics  of  the  energy  containing  eddies.  It  should  be  noted 


that  some  definitions  of  mixing  layer  width  such  as  the  maximum-slope  (or 

vorticity)  thickness  5»(U0-  U,)/(dU/dY)  used  by  Brown  and  Roshko 

w  z  l  max 

(1974),  are  not  sensitive  to  these  variations  in  the  shape  of  the  mean  vel¬ 


ocity  profiles. 


atvv c  ^  w A. »\  •*_ _  . 
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5.2.  Streaklines  and  vorticity  distribution 

One  of  the  main  objectives  of  this  work  was  to  compare  streaklines 
computed  from  phase-locked  components  of  measured  velocity,  with  conven¬ 
tional  flow  vizuallzation  techniques.  Such  a  comparison  could  provide  a 
link  between  qualitative  observations  and  quantitative  data  and  may  serve 
as  a  tool  for  analyzing  the  dynamic  significance  of  the  pairing  of  vor¬ 
tices,  observed  visually  by  Winant  and  Browand  (1974). 

A  comparison  between  a  streakline  pattern  and  a  smoke  picture  is  pre¬ 
sented  in  Fig.  5.2.1  ,  and  the  resemblance  between  the  two  is  obvious.  In 
order  to  take  the  photograph,  a  smoke-filament  was  introduced  near  the  tra¬ 
iling  edge  of  the  splitter-plate,  while  the  test-section  was  illuminated  by 
a  stroboscope  synchronized  with  the  frequency  of  the  flap.  The  exposure 
time  was  0.8  seconds  and,  therefore,  the  photograph  represents  a 
phase-locked  superposition  of  36  events.  The  data  used  for  calculating 
streaklines  was  also  used  for  calculating  the  velocity  and  vorticity 
fields,  the  distribution  of  Reynolds  stresses,  the  transfer  of  energy  from 
the  mean  motion  and  the  scales  of  the  large  coherent  structures.  The  cor¬ 
responding  phase-locked  vorticity  field  is  also  presented  in  this  figure. 

A  pattern  of  4  streaklines,  each  originating  from  a  different 
source-point,  vertically  displaced  at  the  first  measuring-station  (X«200 
mm)  is  shown  at  the  top  of  Fig.  5.2.2  ,  while  the  corresponding 
phase-locked  vorticity  field,  described  by  isodynes  ( iso-vorticity  con¬ 
tours),  is  presented  at  the  bottom.  Each  lumpy  structure  coincides  with  a 


concentration  of  vorticity  and,  therefore,  is  loosely  referred  to  as  a  vor¬ 
tex. 

The  three  regions  discussed  in  conjunction  with  the  spreading  rates  of 
the  shear  layer  can  also  be  distinguished  by  the  manner  in  which  vorticity 
is  distributed.  Vorticity  in  region  I  is  concentrated  in  two  cores,  dis¬ 
placed  both  longitudinally  and  laterally  (Fig.  5.2.3(b)).  The  longitudi¬ 
nal  displacement  of  these  cores  diminishes  with  increasing  distance  down¬ 
stream,  until  it  vanishes  at  the  beginning  of  region  II,  forming  a  circular 
lump.  This  description  agrees  with  predictions  based  on  * nviscid  linear 
stability  theory  and  will  be  discussed  later  in  this  work. 

Region  I,  which  is  marked  by  a  rapid  growth  of  the  shear  layer,  can  be 
characterized  by  the  amplification  of  disturbances  in  the  lateral  direc¬ 
tion.  Its  length  is  equal  to  the  distance  required  for  the  vortex  sheet  to 
roll  into  a  lump.  A  lump  is  formed  at  a  cross  section  at  which  the  streak¬ 
line  pattern  folds  over  backwards.  Fig.  5.2.3(a)  describes  streaklines  in 
region  I. 

The  lateral  dimension  of  the  vortex  lumps  in  region  II  remains  un¬ 
changed  (Fig.  5.2.4(a))  ,  but  their  inclination  to  the  flow  axis  is  al¬ 
tered;  it  crosses  the  vertical  inclination  angle  t/2  at  the  beginning  of 
region  II  (Fig.  5.2.3(a)).  The  thin  vortex  sheet  which  separates  adaja- 
cent  lumps,  and  sometimes  is  referred  to  as  the  "braid",  is  stretched  in 
the  process,  as  can  be  seen  from  the  calculated  streaklines.  The  spatial 
resolution  of  the  measurements  may  be  insufficient  to  detect  the  vorticity 


associated  with  the  braid.  Each  vortex  in  region  II  appears  to  contain  a 
single  core  (Fig.  5.2.4(b)),  as  predicted  by  linear  stability  analysis. 

The  vortex  lumps  in  region  III  are  leaning  downstream  on  the  high  vel¬ 
ocity  side  of  the  shear  layer  (upper-stream  side  of  Fig.  5.2.5(a)),  their 
inclination  angle  to  the  axis  of  the  flow  continuously  decreases  with  in¬ 
creasing  distance  from  the  splitter  plate.  The  elongation  of  the  vortex 
lumps  and  the  change  in  their  inclination  angle  causes  one  lump  to  overlap 
another.  It  seems  that  a  second  vorticity  core  reappears  in  the  interior 
of  each  vortex  (Fig.  5.2.5(b))  in  region  III. 

Neither  the  phase-locked  streaklines  nor  the  photographs  produced  by 
stroboscopic  illumination  show  tha  existance  of  pairing.  Therefore  one 
must  conclude  that  pairing,  if  it  exists,  is  not  locked  in  phase  to  the  ex¬ 
citation  fluctuations  and  it  does  not  occur  very  often  in  a  forced  mixing 
layer. 

Fig.  5.2.6  shows  a  sequence  of  streaklines  plotted  at  time  intervals 
equivalent  to  1/4  period  of  the  forced  oscillations  between  them.  One  can 
therefore  observe  the  temporal  development  in  the  inclination  angle  of  the 
structures  while  moving  from  region  I  to  region  II  and  the  stretching  of 
the  structures  in  region  III.  The  corresponding  vorticity  field  is  pre¬ 
sented  in  Fig.  5.2.7.  The  concentration  ofvorticity  phase-locked  to  the 
frequency  of  excitation,  decreases  with  increasing  X. 
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5.3.  The  intensity  of  the  fluctuations,  the  Reynolds  stress 
and  their  respective  spectral  distributions 

The  inclination  angle  of  the  large  vortex  lump  to  the  direction  of  the 
free  stream  is  related  to  the  sign  of  the  Reynolds  stress  as  schematically 
sketched  below  (Browand  1980).  When  the  inclination  angle  of  the  eddy  Jp is 
greater  than  ft/2  the  associated  Reynolds  stress  is  positive;  when  it  is 
less  than  ft/2  the  Reynolds  stress  changes  sign  and  becomes  negative. 


The  phase-locked  vortices  in  region  I  have  an  elliptical  shape  with  an 
inclination  angle  greater  than  X/2.  It  is  expected,  therefore,  that  the 
phase-locked  -<u,'v,>  product  in  this  region  would  be  positive.  The  oppo¬ 
site  is  true  for  regions  II  and  III  where  the  inclination  angle  of  the  vor¬ 
tices  to  the  mean  flow  is  less  than  X/2  and  -<u7v7>  is  expected  to  be  nega¬ 
tive  there.  The  change  of  sign  occurs  in  the  transition  between  regions  I 
and  II. 

The  Reynolds  st ress  ▼  m-f u*vT  is  linked  to  the  mean  momentum  equation 
and  to  the  mean  kinetic  energy  equation;  therefore,  changes  in  the  distri¬ 
bution  of  -uV'  are  associated  with  changes  in  the  momentum  thickness  ©  and 
with  the  flow  of  energy  to  or  from  the  turbulent  motion. The  overall  energy 
production  term  can  be  subdivided  into  a  phase-locked  (coherent)  term  and  a 


% 
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random  (turbulent)  term.  The  reader  is  referred  to  Hussain's  article 
(1983)  discussing  the  governing  equations  of  the  flow,  using  what  is  gener¬ 
ally  referrred  to  as  the  triple  decomposition. 


The  production  of  energy  integrated  across  the  flow  is  shown  in  Fig. 
5.3.1  as  a  function  of  the  streamwise  distance.  The  overall  production 


term  (marked  by  crosses)  and  the  phase-locked  production  term  (obtained  by 


-<u7v7>^y  and  marked  by  triangles)  are  plotted  together  because  one  may  re¬ 


late  the  energy  production  and  the  Reynolds  stress  to  the  rate  of  spread  of 
the  shear  layer. 


5.3.1.  Region  I 


Transverse  distributions  of  Reynolds  stress  measured  at  four  X  loca¬ 
tions  in  the  portion  of  region  I,  in  which  de/dx  is  constant  (i.e.  X<500 
mm),  are  shown  in  Fig.  5.3.2).  The  data  is  compared  with  the  phase-locked 
-<u’v’>  product,  marked  by  triangles  on  the  same  figure,  which  accounts  for 
most  of  the  Reynolds  stress.  Temporally  averaged  and  phase-locked  energy 
production  terms  behave  in  a  similar  manner  (Fig.  5.3.3),  suggesting  that 
the  flow  is  governed  by  the  phase-locked  fluctuations.  Most  of  the  energy 
(about  75%)  is  transferred  from  the  mean  motion  to  the  phase-locked  fluctu¬ 
ations  that  dominate  the  flow  in  this  region.  Graphs  representing  the 
cross-spectra  at  identical  X  locations  (Fig.  5.3.4),  show  that  the  distur¬ 
bance  frequency  dominates  the  flow,  in  spite  of  the  fact  that  the  first 
harmonic  may  also  be  detected.  Higher  harmonics  do  appear  in  the  power 
spectra  of  the  longitudinal  and  lateral  velocity  components  (Figs.  5.3.5 
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and  5.3.6,  respectively)  but  they  are  not  sufficiently  corrlated  to  influ¬ 
ence  the  Reynolds  stress. 

The  harmonic  distortion,  defined  as  the  ratio  between  the  spectral 
peak  of  the  harmonic  frequency  to  the  excitation  frequency,  may  be  deduced 
from  figures  5.3.5  and  5.3.6.  The  harmonic  distortion  of  the  streamwise 
component  of  velocity  increased  by  approximately  an  order  of  magnitude 
between  X“200  mm  and  X*400  mm  at  the  outer  peripheries  of  the  shear  layer 
but  it  did  not  increase  at  (Y-Yq)/G«0;  the  harmonic  distortion  of  the  la¬ 
teral  velocity  component  increased  by  a  factor  of  2  over  the  same  distance 
and  did  so,  uniformly,  across  the  shear  layer.  The  increase  in  the  harmon¬ 
ic  distortion  with  downstream  distance  suggests  the  ever  increasing  impor¬ 
tance  of  the  non  linear  terms  in  the  equatuion  of  motion  with  increasing  X. 

5.3.2  The  first  transition  subregion 

The  first  transition  subregion  (500<X<640  mm)  is  defined  as  the  region 
in  which  de/dx  changes  from  its  constant  positive  value  in  region  I  to  a 
constant  negative  value  in  region  II.  A  local  maximum  of  the  momentum 
thickness  is  attained  at  X^bO  mm,  where  the  overall  turbulent  energy  pro¬ 
duction  vanishes  (Fig.  5.3.1)  while  the  end  of  this  transition  subregion, 
is  marked  by  a  minimum  of  the  overall  production  (X-64U  mm). 

The  distributions  of  -u'v'  and  <-u'v'>  as  well  as  the  corresponding 
energy  production  terms  resemble  each  other  in  this  subregion  (figs.  5.3.7 
and  5.3.8)  as  might  have  been  expected  from  Fig.  5.3.1.  The  flow  contin- 


ues,  therefore,  to  be  dominated  by  the  excitation  frequency.  The  Reynolds 
stresses  and  the  concomitant  production  terms  change  sign  to  become  nega¬ 
tive.  The  change  is  initially  observed  on  the  high  velocity  side  of  the 
flow  and  is  associated  with  the  increased  inclination  of  the  large  vortices 
resulting  from  the  comparatively  larger  phase  velocity  in  this  region  (see 
Fig.  5-4  in  Gaster,  Kit  &  Wygnanski). 

The  first  harmonic  of  the  cross  spectrum,  which  has  a  magnitude  com¬ 
parable  to  the  amplitude  of  the  imposed  frequency,  changes  its  sign  first 
and  the  imposed  frequency  follows  (Fig.  5.3.9).  The  interaction  between 
the  imposed  fluctuations  and  their  first  harmonic  frequency  appears  to  be  a 
dominant  factor  responsible  for  the  contraction  of  the  mixing  layer  with 
downstream  distance  and  will  be  discussed  later.  (Fig  5.3.9). 

The  power  spectra  of  the  two  velocity  components  in  the  first  transi¬ 
tion  region  are  shown  in  Figs.  5.3.10  and  5.3.11.  The  harmonic  distortion 
of  both  velocity  components  diminishes  in  this  region  particulary  on  the 
high  velocity  side  of  the  mixing  layer.  It  is  Interesting  to  note  that  the 
streamwise  component  of  velocity  at  the  excitation  frequency  is  amplified 
between  X*520  mm  and  X“640  mm  while  the  same  component  at  the  harmonic  fre¬ 
quency  decays;  the  opposite  trend  is  noticable  for  the  transverse  velocity 
component. 

5.3.3  Region  II 

In  region  II  (640<X<840  mm)  the  momentum  thickness  of  the  mixing  layer 


decreses.  The  overall  energy  production  is  negative  (Fig.  5.3.1)  and  the 
Inclination  angle  of  the  vortices  on  the  high-speed-side  of  the  flow  decre¬ 
ases.  The  phase-locked  oscillations  still  dominate  the  flow  (as  can  be 
seen  in  Figs.  5.3.12  and  5.3.13)  and  the  phase-locked  energy  production 
throughout  the  region  is  negative  (Fig.  5.3.1). 

The  fluctuating  velocity  components  in  this  region  are  dominated  by 
the  large  vortices  which  are  perfectly  locked  to  the  phase  of  the  flap  mo¬ 
tion  (see  also  Oster  &  Wygnanski).  Although  the  product  of  <-u*v*>ijY  is 
negative  (Fig.  5.3.13)  suggesting  that  energy  might  be  lost  by  these  large 
eddies  to  the  mean  motion  the  actual  energy  contained  at  the  frequency  of 
excitation  does  not  diminish  between  660<X<840  mm  (Fig.  5.3.15)  the  har¬ 
monic  component  however  (whose  frequency  is  2f)  loses  most  of  its  energy 
between  440<X<840  mm.  The  harmonic  distortion  which  is  0.5  at  X-440  mm  and 
(¥“Y0)/o-2  (Fig.  5.3.5)  is  reduced  to  0.05  at  a  corresponding  lateral  lo¬ 
cation  at  X-840  mm  (Fig.  5.3.15).  The  second  harmonic  (F«3f)  which  is 
significant  in  region  1  disappears  completely  in  region  II  and  most  of  the 

random  fluctuations  are  shifted  towards  frequencies  which  are  lower  than 

2 

the  excitation  frequency.  At  X“200  mm  the  background  spectra  of  u'  in  the 

2 

range  between  0-200  Hz  is  flat  while  the  spectrum  of  v  slopes  upwards  to- 

2 

wards  higher  frequencies.  At  X»440  mm  this  spectrum  of  v'  has  a  broadband 
maximum  around  the  harmonic  frequency  2f  (Fig.  5.3  6);  at  X-640  mm  this 
maximum  shifts  towards  the  excitation  frequency  and  remains  there  till  the 
end  of  region  II  (Fig.  5.3.16).  The  u'v'  spectrum  has  two  negative  peaks 
at  the  excitation  frequency  and  at  its  first  harmonic  between  600<X<800  mm. 
The  relative  importance  of  the  -u^v7  product  at  2f  decreases  with  increas- 
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ing  X  in  this  region  until  it  vanishes  across  the  entire  mixing  layer  at 
X-840  mm  (Figs.  5.3.9  and  5.3.14);  the  <-u,vT>  product  at  the  excitation 
frequency  becomes  negative  at  X-600  mm  and  remains  negative  until  X**900 
mm). 


5.3.4  The  second  transition  region 

The  transition  from  region  II  to  III  (840<X<1080  mm)  is  marked  by  a 
renewed  growth  of  the  momentum  thickness  with  Increasing  X  and  a  positive 
overall  production  of  turbulent  energy. 

In  this  region  we  observe,  for  the  first  time,  that  the  phase-locked 
-<u’v’>  product  and  the  associated  production  of  energy  deviate  from  the 
time-averaged  -u’v’  product  and  overall  energy  production  (Figs.  5.3.17 
and  5.3.18);  while  the  overall  production  is  positive,  the  phase-locked 
production  of  energy  is  negligibly  small.  Positive  values  do  appear  in  the 
phase-locked  -<u,’v’>  product  and  <-u’v,>j[y  in  the  centeral  region  of  the 
mixing  layer,  but  their  influence  is  offset  by  the  negative  values  at  the 
peripheries.  Positive  peaks  in  the  cross  spectra  appear  in  both  the  dis¬ 
turbance  frequency  and  its  first  harmonic  in  the  vicinity  of  (Y-Yq)/$-o. 
(Fig  5.3.19). 

2  2 

The  spectral  distributions  of  u  and  v  do  not  change  appreciably 

between  880<X<1080  mm;  the  harmonic  distortion  coefficient  remains  con- 

2 

stant  and  the  background  turbulence  of  the  v'  fluctuations  has  a  broadband 
maximum  around  the  excitation  frequency  (Figs  5.3.20  amd  5.3.21). 
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Since  the  momentum  thickness  in  this  region  is  smaller  than  the  maxi¬ 
mum  6  attained  at  the  end  of  region  I,  the  mixing  layer  may  return  to  be 
linearly  unstable  to  the  excitation  frequency.  This  may  be  the  reson  for 
the  reappearance  of  positive  phase-locked  -<u*v*>  and  the  corresponding 
production  terms. 

5.3.5  Region  III 

In  region  III  (X>1080  mm)  the  mixing  layer  resumes  its  linear  rate  of 
spread  with  X.  The  growth  rate,  however,  is  approximately  equal  one  half 
of  the  growth  rate  in  region  I.  The  overall  energy  production  in  this  re¬ 
gion  is  positive  everywhere  (Fig.  5.3.1),  and  both  -u'v'  product  and  over¬ 
all  energy  production  profiles  are  positive  (Figs.  5.3.22  and  5.3.23,  res¬ 
pectively).  On  the  other  hand,  the  phase-locked  -<u’vT5  and  -<u*v,>^-  dis¬ 
tributions  are  negative  in  this  region,  consequently  (-<u*vT>4^-  dY  is  ne- 
gative  (Fig.  5.3.1). 

The  growth  of  the  mixing  layer  does  not  appear  to  be  controlled  by  the 
external  excitation,  since  the  phase-locked  vortices  are  losing  energy 
through  the  “<uyv,>-gy  product.  The  observation  is  reinforced  by  the  nega¬ 
tive  cross  spectral  peaks  occuring  at  the  excitation  frequency  (Fig 
5.3.24),  however,  a  broad  band  of  positive  -u*v*  product  appears  at  fre¬ 
quencies  which  are  lower  than  the  imposed  frequency.  In  fact,  the  intensi¬ 
ty  of  the  fluctuations  at  the  excitation  frequency  decreases  very  slowly 
with  downstream  distance  for  X>1100  mm  but  the  background  low  frequencies 
increase  by  a  factor  of  5  between  1100<X<1700  mm  (Figs.  5.3.25  and 
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5.3.26).  Therefore,  eddies  which  are  larger  than  the  eddies  associated 
with  the  excitation  frequency  dominate  this  region.  These  eddies  are  not 
locked  in  phase  to  the  imposed  oscillations  and  are  presumably  responsible 
for  the  resumed  growth  of  the  shear-layer. 

5.3.6  The  subha rmonic  frequency 

Anticipating  that  the  resumed  growth  of  the  mixing  layer  in  region  III 
stems  from  vortex  pairing  special  attention  was  paid  to  this  process  and  to 
the  generation  of  the  subharmonic  frequencies.  Since  the  large  eddies 
passing  through  region  II  were  locked  to  the  phase  of  the  excitation  it  was 
assumed  that  the  pairing  process  will  follow  suit  and  therefore,  the  phase 
averaging  was  done  at  the  subharmonic  frequency  as  described  in  chapter  4. 
A  comparison  between  the  data  phase  locked  to  the  excitation  signal  or  its 
subharmonic  could  not  be  distinguished  suggesting  that  the  process,  if  it 
exists,  occurs  randomly  in  time  or  space.  Power  spectra  of  u  and  v  re¬ 
veal  interesting  exchanges  between  the  excitation  frequency  and  its  harmon¬ 
ic  component  (F"2f)  in  regions  I  and  II  and  a  good  portion  of  region  III 
(up  to  X»1300  mm)  but  they  do  not  show  a  clear  appearance  of  a  subharmonic 
frequency  (F-f/2)  as  was  initially  expected  (Fig  5.3.26).  Comparing,  say, 
the  last  spectrum  shown  in  Fig.  5.3.1b  (taken  at  X-840  mm  and  (Y-Yq)/©*2) 
to  a  corresponding  spectrum  in  Fig.  5.3.26  (X-1700  mm)  one  notices  the 
following  (i)  disappearance  of  peaks  at  the  harmonic  frequency  at  X*1700  mm 
(ii)  the  disappearance  of  side  bands  (i.e.  the  broadening  of  the  spectrum) 
at  the  excitation  frequency  (ill)  the  decay  of  the  peak  at  the  peak  at  the 
excitation  frequency  (iv)  the  increase  in  the  background  level  of  all  fre- 
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quencies  that  F<£  by  approaximatly  an  order  of  magnitude. 


Cross  spectra  shown  in  Fig.  5.3.24,  smoke  pictures  illuminated  at  the 
subharmonic  frequency,  phase-locked  streakline  patterns  and  vorticity  con¬ 
tours  did  not  show  any  noticeable  activity  associated  with  pairing.  It  was 
therefore  concluded  that  high  amplitude  excitation  might  have  inhibited  the 
pairing  process,  as  we  understood  it.  We  are,  therefore,  continuing  with 
experiments  in  which  the  flow  is  excited  by  two  frequencies  which  are 
phase-locked  to  one  another  hoping  to  lock  onto  the  pairing  process. 


* 

I 


CHAPTER  6 


A  DETAILED  COMPARISON  WITH  THE  LINEAR  INVISCID  STABILITY  THEORY 


The  inviscid  stability  theory  provides  a  powerful  tool  In  modeling  the 
development  of  slightly  disturbed  flows.  Gaster,  Kit  and  Wygnanski  (to  be 
published  in  J.F.M.)  applied  this  model  to  a  turbulent,  slowly  diverging, 
mixing  layer  and  compared  the  theoretical  calculations  for  the  streamwise 
velocity  component  with  experimental  results.  Although  the  lateral  distri¬ 
butions  of  normalized  amplitudes  and  phases  showed  a  good  agreement  with 
the  model,  the  predicted  amplification  rate  exceeded  the  experimental  ob¬ 
servations.  A  partial  purpose  of  this  experiment  was  to  increase  the  exci¬ 
tation  level  until  a  harmonic  distortion  occurs  and  observe  the  shortcmo- 
nigs  of  the  linear  theory.  The  measurements  were  also  extended  to  Include 
two  velocity  components  from  which  phase-locked  vorticity  and  -u’v’  product 
were  calculated.  All  measurements  were  compared  with  the  theoretical  cal¬ 
culations. 

The  inviscid  equation  of  motion  describing  the  vorticity  of  a 
two-dimensional  flow  can  be  written  as  follows: 

(6.1) 

Et  at  “a*  3n 

2-  2 

where  n_V  f+V  T  ; 

y( £ , n )  is  a  complex  two-dimensional  mean  stream  function; 

is  a  complex  time-dependent  stream  function; 

£  and  ^  being  the  dimensionless  longitudinal  and  lateral  coordinates,  res- 


where  the  slow  variation  with  £  is  taken  up  by  the  amplitude  A(()  and  the 
eigenfunction  «(£,n)  (i.e.  their  first  derivatives  with  respect  to  £  are 

of  order  e)  •  Substituting  (6.2)  into  (6.1)  and  omitting  higher  order 
terms  (i.e.  second  derivatives  of  *  and  A  with  respect  to  (  or 
cross-products  of  first  order  derivatives),  one  obtains  after  some  manipu¬ 
lation  two  equations:  (i)  the  local  inviscid  Orr-S ommerf eld  equation  for 
the  equivalent  parallel  flow,  involving  the  local  mean  velocity  profile  for 
order  1: 

(6.3)  iA{(*"-a2*)(Ua-w)-U"*a}-0 

and  (ii)  an  equation  for  order  e,  from  which  the  amplitude  A(fc)  can  be  cal¬ 
culated: 

A(e)-exp[-J|^  dx] 

where  N(()  and  M(()  are  the  following  integrals: 

m 

2aw*+U[*"-3a2*]-U"*}*avl 

N(C)“  /U[*(3a/3C)  +  2a(3$/H)]  +  U[  (  3*"/h)-302(  3*/3e)~3a*(  3«/3  O ) 

+  *'(3U'/3t)-U"(3*/35)  +  V[*"'-a2*']}^ 


The  notation  ,  represents  the  first,  second  and  third  deriva¬ 
tives  with  respect  to  n ,  respectively.  This  equation  stems  from  a  condi¬ 
tion  described  by  Ince  (1956):  "A  necessary  and  sufficient  condition  that 
a  non-homogeneous  linear  system  with  homogeneous  boundary  conditions  may 
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have  a  solution  is  that  every  solution  of  the  homogeneous  adjoint  system 
(homogeneous  equation  and  boundary  conditions)  satisfies  the  relation:" 

•O 

Jg$dn«0 

-o» 

This  condition  is  commonly  referred  to  as  the  solvability  condition.  The 
function  ♦  is  the  eigenfunction  corresponding  to  an  eigenvalue  of  the  adjo¬ 
int  Orr-S ommerfeld  system  and  has  the  form: 

^“♦/(0-w/a) 

The  shape  of  the  mean  velocity  profile  used  for  the  solution  of  the 
problem  was  a  fit  based  upon  the  hyperbolic  tangent  profile  and  its  deriva¬ 
tives  (suggested  by  Prof.  R.  E.  Kaplan  from  the  University  of  Southern 
California): 

U(n/?)  •  0.5  U+Tanh(n/£)  [l-C/cosh2  (n/OJ} 

The  constant  C  was  chosen  to  be  -0.67  .  This  profile  fits  the  experimental 
data  better  than  the  hyperbolic  tangent  profile  commonly  used  in  the  liter¬ 
ature.  The  chosen  analytical  expression  is  marked  by  dots  in  Fig.  6.1, 
while  the  experimental  profiles  are  marked  by  a  variety  of  symbols. 

The  method  for  calculating  the  theoretical  eigenfunctions  is  described 
in  detail  by  Crighton  &  Gaster  (1976)  and  by  Gaster,  Kit  &  Wygnanski.  The 
dimensionless  coordinates  for  the  experimental  results  were  calculated  by: 

t  -  <x-*0>a!  /  «t  /  i 

<*-v  1 8 

where  6t  represents  the  theoretical  momentum  thickness  at  a^»l  (9t*.2467); 

Xq  is  the  virtual  origin  of  the  flow  where  the  momentum  thickness  vanishes 

da 

•JU  yields  the  growth- rate  of  the  experimental  momentum  thickness; 


Y0  is  the  lateral  coordinate  corresponding  to  the  location  at  which 
0»O.5(Uj+U2); 

6  is  the  local  experimental  momentum  thickness. 

The  definition  of  L  is:  L-(U2"U1  )CQ/2irf  ; 
where  Cq  is  a  constant  chosen  to  be  3  in  order  to  ensure  that  a^-o  at  the 
end  of  the  calculation  domain; 

f  is  the  experimental  excitation  frequency  in  Hz 

(U2~U^)  is  the  difference  between  the  mean  free-stream  velocities  in 
mm/sec. 

The  two  components  of  the  oscillatory  motion  described  by  amplitude 

and  phase  were  calculated  by  Fourier  representation  of  the  velocity  signals 

locked  to  the  disturbence  frequency.  Two  coefficients  were  calculated  for 

each  phase-locked  signal: 

x  r 

A*^J(u(t)8in(u)yt)dt  and  B-4^u(t)cos(w0t)dt 

O  0 

(where  Wy  is  the  frequency  of  the  fluctuations),  or  using  the  numerical 
scheme: 

NPP  NPP 

A-  Z  <u( j)>sin[(j-l)*2*/NPP)J  and  B-  <u( j)>cos[ ( j-l)*2*/NPP) ] 

J-l 

The  amplitude  and  phase  of  each  velocity  signal  were  derived  by: 

r 2  r 

O  y  A  +B  and  $"arctan(A/B) 

The  -u’v’  product  amplitude  is  by  definition: 

T“C1*C2*CC6(W 

where  the  subscripts  1  and  2  stand  for  the  longitudinal  and  lateral  direc¬ 
tions,  respectively. 
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The  amplitude  and  the  phase  distributions  of  the  spanwise  vorticity 
component,  were  computed  by  numerical  differentiation  of  neighboring 
phase-locked  velocity  vectors,  which  were  then  Fourier  transformed. 

The  comparison  between  the  experimental  results  and  the  theoretical 
calculations  is  divided  into  two  categories: 

A  comparison  between  normalized  lateral  distributions  of  the  following 
functions: 

a)  The  amplitude  of  the  longitudinal  velocity  component 
(Figs.  6. 2(a), (b)). 

b)  The  phase  of  the  longitudinal  velocity  component 
(Figs.  6. 3(a), (b)). 

c)  The  amplitude  of  the  transverse  velocity  component 
(Figs.  6. 4(a), (b)). 

d)  The  phase  of  the  transverse  velocity  component 
(Figs.  6. 5(a), (b)). 

e)  The  amplitude  of  the  vorticity  (Figs.  6. 6(a), (b)) 

f)  The  phase  of  the  vorticity  (Figs.  6. 7(a), (b)) 

g)  The  amplitude  of  -u^  product  (Figs.  6. 8(a), (b)) 

The  triangles  in  all  the  figures  represent  experimental  results  while 
the  theoretical  predictions  at  corresponding  £  stations  are  marked  by  solid 
curves.  Both  the  dimensionless  coordinate  £  and  the  corresponding  stream- 
wise  distance  from  the  virtual  origin  appear  above  each  curve. 
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•  The  amplitude  distribution  of  any  quantity  Q  in  Figs.  6.2-6.#  was 

m 

normalized  by  the  local  integral  ^|Q|dn>  The  phase  distributions  were 
matched  at  a  single  arbitrary  point ,  namely  at  n“0. 

\ 

’  The  good  agreement  between  all  the  calculated  and  measured  quantities 

f  indicates  that  the  linear  inviscid  theory  serves  as  a  good  first  step  in 

predicting  the  growth  of  the  disturbances  in  region  I  of  a  turbulent  mixing 
layer. 


The  overall  amplification  of  the  longitudinal  velocity  in  the  stream- 
wise  direction  which  is  represented  by  J |u£|dn  is  shown  in  Fig.  6.9, the 

-  mo 

corresponding  integral  for  lateral  velocity  component  is  plotted  in  Fig. 

6.1 0  while  the  spanwlse  vorticity  and  |-u'v'|  product  are  shown  in  Figs. 

6.11  and  6.12,  respectively.  All  the  figures  are  plotted  on 

semi-logarithmic  coordinates  in  which  the  abscissa  is  the  dimensionless 
streamwise  coordinate  and  the  ordinate  is  the  local  value  of  the  integral 
normalized  by  its  value  at  the  first  station.  The  triangles  represent  the 
experimental  data  while  the  theoretical  values  are  plotted  by  solid  curves. 


All  four  integrals  calculated  theoretically  increase  with  X,  the  meas¬ 
ured  amplitudes  of  the  two  velocity  components  also  increase  with  X,  but  at 
considerably  slower  rate.  The  experimental  and  the  theoretical  integrals 
of  the  fluctuations  in  the  lateral  direction  increase  much  more  rapidly 
with  X  than  the  streamwise  component  of  the  fluctuations  (cf.  Fig.  6.9 
with  Fig.  6.1U).  The  integarl  of  the  experimental  longitudinal  velocity 
component  tends  to  decrease  towards  the  end  of  the  domain  of  comparison. 


The  vorticity  amplitude,  integrated  across  the  flow,  decreases  with  X,  and 
the  integrated  -u*v*  product  remains  constant  over  most  of  the  domain  of 
comparison  and  decreases  towards  the  end  of  this  domain.  The  discrepancy 
between  the  high  amplification  rates  predicted  by  the  linear  model  and  the 
relatively  low  amplification  rates  realized  in  the  experiment  deserves  spe¬ 
cial  attention  in  view  of  the  excellent  agreement  between  the  shapes  of  the 
theoretical  and  experimental  eigenfunctions. 

The  linear  model  used  assumes  that  the  flow  departs  a  little  from  the 
parallel  flow  geometry,  and  a  weak  modification  in  the  parallel  stabilty 
theory  is  sufficient  to  describe  the  resulting  solution.  The  small  scale 
turbulence  was  not  accounted  for  by  the  linear  model,  but  its  existance  may 
contribute  somewhat  to  the  above-mentioned  discrepancy  between  theory  and 
experiment.  The  turbulent  energy  production  in  region  I  accounts  for  25% 
of  the  total  (Fig.  5.3.1),  while  the  remaining  75%  is  attributed  to  the 
phase-locked  fluctuations  at  the  disturbence  frequency  and  its  first  har¬ 
monic.  The  contribution  of  the  first  harmonic  was  not  accounted  for  in  the 
analysis  of  the  data  but  it  is  not  insignificant  (Fig.  6,13) 

Another  aspect  that  must  be  considered  when  applying  a  linear  model  is 
the  level  of  the  fluctuation.  The  displacement  amplitude  of  the  oscillat¬ 
ing  flap  in  the  experiment  described  by  Gaster,  Kit  &  Wygnanski  was  much 
smaller  than  the  amplitude  in  the  present  experiment.  In  spite  of  it,  they 
found  that  45%  of  the  overall  energy  at  £“0.786  occured  at  the  imposed  fre¬ 
quency,  while  the  maximum  level  of  the  streamwise  fluctuations  at  this 
point  was  0.18.  These  fluctuation  levels  are  not  small  and  it  is  not  at 
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all  obvious  that  the  linearization  of  the  equations  of  motion  is  justified. 
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CHAPTER  7 

CONCLUDING  REMARKS 


The  strongly  excited  turbulent  mixing  layer  was  examined  In  great  de¬ 
tail,  nevertheless,  many  questions  remain  to  be  answered.  For  example: 

(1)  Can  the  big  Initial  growth-rate  of  the  6hear  layer  be  predicted? 

(2)  What  causes  the  momentum  thickness  to  level  off  at  the  end  of  region 
I  and  what  is  the  cause  for  the  resumption  of  growth  in  region  Ill? 

(3)  Why  does  the  shear  layer  actually  shrink  in  region  II  when  the  level 
of  excitation  is  high? 

(4)  Why  do  the  predicted  shapes  of  the  eigenfunctions  agree  so  well  with 
experiment  and  the  rates  of  amplification  do  not? 

We  believe  that  some  of  the  answers  can  be  provided  by  expanding  the 
linear  calculations  to  account  for  the  nonlinear  processes.  The  perturba¬ 
tions  introduced  by  the  flap  have  a  sizeable  harmonic  component  (at  a  fre¬ 
quency  F**2f)which  is  initially  amplified  with  X  more  swifly  than  the  exci¬ 
tation  frequency  itself.  As  the  mean  flow  broadens  and  becomes  neutrally 
stable  to  the  harmonic  frequency,  on  the  basis  of  the  inviscid  linear 
model,  the  amplitude  of  the  harmonic  fluctuations  diminishes.  The  process 
starts  to  occur  around  X«500  mm  in  this  case,  where  the  -u*v"'  product  of 
F*2f  becomes  negative  (Fig.  5.3.9)  nevertheless  remanents  of  the  harmonic 
frequency  are  still  noticeable  at  X“130Q  mm.  The  slow  decay  of  the  harmon¬ 
ic  amplitude  may  be  attributed  to  a  contribution  from  the  fundamental  ire- 


quency  through  a  nonlinear  process. The  exchange  of  energy  between  the  two 
frequencies  is  currently  investigated  in  a  fully  controlled  experiment. 

The  -u*v*  product  at  the  excitation  frequency  changes  sign  at  X*560  mm 
corresponding  to  fe/Uc-0.030  which  is  the  neutrally  stable  Strouhal  number 
for  the  local  mean  velocity  profile.  The  subharmonic  frequency  at  this  lo¬ 
cation  was  insignificant  and  could  not  sustain  the  amplitude  of  the  of  the 
fundamental  frequency  by  nonlinear  means;  some  exchange  could  have  taken 
place,  however,  between  the  harmonic,  the  fundamental  and  the  mean.  The 
negative  -u7v7  product  at  the  fundamental  frequency  (Fig.  5.3.9)  affected 
the  mean  flow  field  causing  an  actual  shrinkage  of  the  shear  layer.  This 

observation  is  consistent  with  the  observations  of  Oster  &  Uygnanski  (1982, 

d.0 

Fig.  12)  who  also  noticed  that  a  high  amplitude  oscillation  produced  <o 
in  region  II.  At  the  end  of  region  II  (at  900<X<1100  mm)  the  shear  layer 
became  sufficiently  narrow  and  the  shape  of  the  mean  velocity  was  suffi¬ 
ciently  altered  (Fig  5.1.3)  that  the  shear  layer  became  again  unstable  to 
the  fundamental  frequency  (Fig.  5.3.19).  At  this  location,  the  amplitudes 
of  the  random  fluctuations  at  low  frequencies  became  much  larger  and  the 
shear  layer  resumed  its  slow  growth. 

The  explanation  put  forward  for  the  growth  of  the  mixing  layer  is 
based  on  the  premise  that  waves  interacting  in  a  nonlinear  manner  produce 
the  observed  behaviour  of  the  excited  shear  layer.  It  may  serve  as  an  al¬ 
ternative  explanation  to  the  pairing  process  provided  it  will  be  proven  ri¬ 
gorously  and  supplemented  by  additional  experimental  verification. 
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Figure 

1.1  Schematic  sketch  of  the  two-dimensional  mixing  layer. 

2.1  The  wind  tunnel. 

5.1.1  The  variation  of  the  momentum  thickness ,  ©  ,  and  the  width  of 
the  mixing  layer,  b,  with  distance  from  the  splitter  plate. 

5.1.2  The  variation  of  YQ  1Q,  YQ  5Q  and  YQ  g5  with  distance  from 

the  splitter  plate.' 

5.1.3  A  comparison  between  the  normalized  mean  velocity  profiles 
in  regions  I  and  III. 

5.2.1  The  phase-locked  structures:  (a)  A  smoke  picture  with 

illuminated  stroboscopically.  Exposure  time  -  0.8  sec.  (b) 
Calculated  streaklines,  (c)  Calculated  isodynes. 

5.2.2  The  phase-locked  structures:  (a)  Calculated  streaklines, 

(b)  Calculated  isodynes. 

5.2.3  The  phase-locked  structures  in  region  I:  (a)  Calculated 

streaklines,  (b)  Calculated  isodynes. 

5.2.4  The  phase-locked  structures  in  region  II:  (a)  Calculated 

streaklines,  (b)  Calculated  isodynes. 

5.2.5  The  phase-locked  structures  in  region  III:  (a)  Calculated 

streaklines,  (b)  Calculated  isodynes. 

5.2.6  Temporal  evolution  of  streaklines. 

5.2.7  Temporal  evolution  of  isodynes. 

5.3.1  The  variation  of  laterally  integrated  overall  and 

phase-locked  energy  production  with  X. 

5.3.2  The  distribution  of  -u*v9/ (U^U^ )2  and  -<u'vT>/(U2-U.)2  in 

region  I. 

5.3.3  The  distribution  of  -^~v?^'V(U9-U1  )3  and 

-<u’vT>^V(U2-U1)-i  in  region  I.  dY  1  1 

5.3.4  The  cross  spectra  of  u'  amd  v'  in  region  I. 

5.3.5  The  power  spectra  of  the  streamwise  velocity  component  in 

region  I. 


5.3.6  The  power  spectra  of  the  lateral  velocity  component  in 
region  I. 
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5.3.7  The  distribution  of  -u'v'/d^-l^)2  and  -<u'v'>/(U2~U1  )2  in 

the  first  transition  region. 

n  of  -u'v 

in  the  first  trans 

5.3.9  The  cross  spectra  of  u'  amd  v'  in  the  first  transition 

region. 

5.3.10  The  power  spectra  of  the  streamwise  velocity  component  in 

the  first  transition  region. 

5.3.11  The  power  spectra  of  the  lateral  velocity  component  in  the 
first  transition  region. 

5.3.12  The  distribution  of  -u'v'/(U2-U. )2  and  -<u'v'>/(U2-U,)2  in 

region  II. 

5.3.13  The  distribution  of  -^v,4H*/(U,-U1  )3  and 

-<S’7’>gV(U2-01)3  in  region  II.  **  1  1 

5.3.14  The  cross  spectra  of  u'  amd  v'  in  region  II. 


5.3.15 

The  power  spectra  of  the  streamwise 
region  II. 

velocity 

component 

in  { 

1 

5.3.16 

The  power  spectra  of  the  lateral 
region  II. 

velocity 

component 

in  J 

5.3.17 

The  distribution  of  -u'v'/ (U  -U, )2  and  -<u'v' 

>/(U  -u.)2 

in  j 

the  second  transition  region. 

4  1 

I 

5.3.18 

The  distribution  of  -u 

^/(u2 

-u.)3 

and 

-<u'v'>jY*/(U2“U|y  in  the  second  transition  region. 

5.3.19  The  cross  spectra  of  u'  amd  v'  in  the  second  transition 
region. 

5.3.20  The  power  spectra  of  the  streamwise  velocity  component  in 

the  second  transition  region. 

5.3.21  The  power  spectra  of  the  lateral  velocity  component  in  the 
second  transition  region. 

5.3.22  The  distribution  of  -u^v*/ (b2-U. )2  and  -<utvt>/(U2~IK  )2  in 

region  III.  1 

5.3.23  The  distribution  of  and 

-<u'v'>^y7/(U2~U^)3  in  region  III. 

5.3.24  The  cross  spectra  of  u'  amd  v'  in  region  III. 

5.3.25  The  power  spectra  of  the  streamwise  velocity  component  in 


5.3.8  The  .distributio 


"dU  v . . 

ition  region. 
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m 

|  region  III. 


j  5.3.26 

The  power  spectra  of  the  lateral  velocity  component  in 
region  III. 

jj  6.1 

The  normalized  mean  velocity  profiles  in  region  I. 

J  6.2 

The  lateral  distribution  of  the  amplitudes  of  u£  at 
streamwise  locations. 

various 

•  6.3 
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.a 

The  lateral  distribution  of  the  phase  angles  of 
various  streamwise  locations. 

u'  at 

1 

•  6-4 

* 

j 

The  lateral  distribution  of  the  amplitudes  of  v'  at 
streamwise  locations. 

various 

!  6.5 

The  lateral  distribution  of  the  phase  angles  of 
various  streamwise  locations. 

v'  at 

;  6.6 

t 

[ 

The  lateral  distribution  of  the  amplitudes  of  the 
vorticity  component  at  the  excitation  frequency  at 
streamwise  locations. 

spanwise 

various 

■:  .  6.7 

- 

• 

- 

■„ 

The  lateral  distribution  of  the  phase  angles  of  the 
vorticity  component  at  various  streamwise  locations. 

spanwise 

\  6.8 

t 

The  lateral  distribution  of  amplitudes  of  -<u’vT> 
excitation  frequency  at  various  streamwise  locations 

at  the 

• 

l  6.9 

V 

V 

The  amplification  of  the  laterally  integrated  longitudinal 
disturbances  with  distance  from  the  splitter  plate. 

J  6.10 

The  amplification  of  the  laterally  integrated 
disturbances  with  distance  from  the  splitter  plate. 

lateral 

t  6. H 

The  amplification  of  the  laterally  integrated  spanwise 
vorticity  component  with  distance  from  the  splitter  plate. 

C*>J 

• 

The  amplification  of  the  laterally  integrated  -u'v' 
with  distance  from  the  splitter  plate. 

product 

J  6.13 

The  variation  of  J|u£|dY  and  J|u^£-|dy  with  downstream 

distance..  "**  -«• 
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Fig.  5.3. 
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Fig.  5.3.5 
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AD-A179  717 


UNCLASSIFIED 


A  HIQHLV  EXCITED  TUABULENT  MIXING  LAVEIt<U>  TEL-AVIV 
UNIV  (ISRAEL)  DEPT  OF  FLUID  MECHANICS  AND  HEAT  TRANSFER 
I  HEISBROT  ET  AL.  NOV  B4  AE-B4-2  AF0SR-TR-87-8344 
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*'  -■ 
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X*28Q  mm  X“30Q  mm 
^•-0.  493  ^-0.  543 


380  mm 
0.  743 


X-200  mm  X-220  mm  X-240  mm  X-260  mm  X-280  mm  X-300  mm 

r-o.  293  X“°- 343  T“°* 393  T“°- 443  T"°- 493  543 


NORMALIZED  AMPLITUDES  OF  V 


X— 320  mm  x=340  mm  X=360  mm  X=380  mm  Xas400  mm  X*420  mm 

^*=0.  593  f=0.  643  T“0-  693  \=0. 743  'po.  793  T=0- 843 


NORMALIZED  AMPLITUDES  OF 


'A  C’ATa 


210  mm  X-230  mm  I  X-250  mm  I  X-270  mm  I  X-290  mm  I  X— 310  mm 

0.318  r-0. 368  f“0.  418  5-0. 468  T-Q.  518  5-0.568 


AMPLITUDES  OF  VQRTICITY 


0.618  T-0.668  f-0. 718  f 0. 768  ^-0.818  f-0. 868 


PHASE  OF  VORTICITY 
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PHASE  OF  VORTICITY 


NORMALIZED  AMPLITUDES 
OF  -IFV7  PRODUCT 

Fig.  6.  8  (a) 


320  mm  X-340  mm  X-360  mm  X-380  mm  X-400  mm  X-420  mm 

0.593  f-0.  643  f  «0. 693  T-0. 743  f*0. 793  f-0. 843 
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OF  -U'V'  PRODUCT 
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